Activated Random Walkers: 
Facts, Conjectures and Challenges 

Ronald Dickman, Leonardo T. RoUa, Vladas Sidoravicius 
December 29, 2009 

Abstract 

We study a particle system with hopping (random walk) dynamics on the integer 
lattice Tf" . The particles can exist in two states, active or inactive (sleeping); only 
the former can hop. The dynamics conserves the number of particles; there is no 
limit on the number of particles at a given site. Isolated active particles fall asleep 
at rate A > 0, and then remain asleep until joined by another particle at the same 
site. The state in which all particles are inactive is absorbing. Whether activity 
continues at long times depends on the relation between the particle density C and 
the sleeping rate A. We discuss the general case, and then, for the one-dimensional 
totally asymmetric case, study the phase transition between an active phase (for 
sufficiently large particle densities and/or small A) and an absorbing one. We 
also present arguments regarding the asymptotic mean hopping velocity in the 
active phase, the rate of fixation in the absorbing phase, and survival of the infinite 
system at criticality. Using mean-field theory and Monte Carlo simulation, we 
locate the phase boundary. The phase transition appears to be continuous in both 
the symmetric and asymmetric versions of the process, but the critical behavior is 
very different. The former case is characterized by simple integer or rational values 
for critical exponents (/3 = 1, for example), and the phase diagram is in accord with 
the prediction of mean-field theory. We present evidence that the symmetric version 
belongs to the universality class of conserved stochastic sandpiles, also known as 
conserved directed percolation. Simulations also reveal an interesting transient 
phenomenon of damped oscillations in the activity density. 

KEYWORDS: Interacting particle systems; absorbing-state phase transition; sandpiles; 
interacting random walkers 

1 Introduction 



Interacting particle systems with conservation have attracted great interest in physics, 
probability, and allied fields, in part because they afford simple examples of phase tran- 
sitions in systems maintained far from equilibrium. In these models the local dynamics 
conserves the number of particles, although certain sites may act as particle sources or 
absorbers. One broad important class of models subsumes exclusion models, in which 
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particles interacting via on-site exclusion (and possibly an additional short-range inter- 
action) execute biased hopping on a lattice. Important examples are driven diffusive 
systems [30l [HI [SU |44] and the totally asymmetric exclusion process [52l [37l H] . In 
another class of models there is no exclusion (any number of particles may occupy the 
same site) but the particles exist in two states that may be termed active and inactive, 
such that activation of an inactive particle requires the intervention of one or more 
active ones. This class includes so-called conserved lattice gases [50l [38l [Ml [iQl [H] 
and stochastic sandpile models [l2l [ISj [lU [El [TO]. Such models exhibit self-organized 
criticality [21 [3l [231 [S] when coupled to a suitable control mechanism |13i (8j . 

In this paper we study a system of activated random walkers (ARW) on the lattice. For 
theoretical analysis, it is convenient to define the model (ARWl) on the infinite integer 
lattice Z^. We assume, in this case, that there are infinitely many particles in the 
system, each of which can be in one of two states: A (active) or S (inactive or sleeping). 
Each 74-particle performs an independent, continuous time, simple symmetric random 
walk on Z*^, with the same jump rate, which we assume, without loss of generality, to be 
equal to 1. When an ^-particle jumps to a site with an S-particle or particles, any such 
particle at this site is immediately activated (i.e., switches to state A). Each isolated 
^-particle goes to sleep (switches to state 5), at a rate A > 0. [From this rule it follows 
that if two or more particles occupy the same site, then they are all of type A or all of 
type S (the latter situation can only arise in the initial condition).] Since 5-particles are 
immobile, at any given site, at most one A-particle can go to sleep, and (if undisturbed), 
remain in state S forever after. The limit A — > oo corresponds to the model studied by 
Jain [26], in which any isolated particle immediately becomes immobile. 

We assume that initially the particles are distributed according to a product Poisson 
measure with mean and all the particles are active. 

In numerical studies the following one-dimensional model (ARW2) is used. The system 
is a chain of L sites with either periodic or open boundaries. Initially N particles are 
randomly placed in the system. (In the case of periodic boundaries the particle number 
is conserved.) As in ARWl, each nonisolated A-particle hops at unit rate. An isolated 
^-particle has a somewhat different dynamics: it hops at rate p < 1 and goes to sleep 
at rate q = 1 — p. Thus in ARW2 isolated ^-particles have a smaller hopping rate 
than nonisolated ones, while in the ARWl all A-particles have the same jump rate. 
While certain details such as the phase boundary may differ between the two versions, 
we expect the global properties to be the same. For the question of whether or not 
the system fixates, the two models are equivalent via A = q/p. The case p = again 
represents the model studied in [26] . 




A natural generalization of the model is to allow the ^-particles to execute a biased 
random walk. We shall in fact be particularly interested in the completely asymmetric 
case. Another possible generalization is to assume that each 5-particle is activated at 
rate < a < +oo, when it shares its site with an A-particle. One could go even further, 
by assuming that the rate at which an A-particle activates S particles at a given site 
depends on the number of S particles at the site (zero-range rule). There is however a 
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substantial difference between the case a < +00, and tlie case a = +00, described at 
tlie beginning of this section. If a < +00, during evolution the ^-particles may share a 
site along with one or more S'-particles, so that this model is a kind of dynamic contact 
process. If a = +00, the situation changes, and we believe this model belongs to the 
universality class of stochastic conserved sandpiles. In this work, we only consider the 
case a = +00 (instantaneous reactivation). 

The primary motivation for the present study is the stochastic conserved sandpile, 
generally known as Manna's model |42i I43j . In infinite volume, this model is defined 
as follows. We assume that initially there are infinitely many particles, distributed in 
such a way that at each site of Z*^ we have Poisson mean C > number of particles. 
Each site is equipped with an exponential rate 1 clock, and each time a clock rings at 
a site bearing 2d or more particles, 2d particles move from the site to randomly chosen 
nearest neighbors. That is, (differently from the deterministic Bak-Tang-Wiesenfeld 
sandpile model [21 [3]), each particle chooses its direction among the 2d possibilities with 
probability {2d)~^ , independent of any other particle. In contrast to the deterministic 
sandpile [6], very little is known rigorously about this system, and the ARW model is a 
reasonable caricature that seems to capture some essential aspects of Manna's model. 

The ARW model may also be viewed as a special case of a diffusive epidemic process, 
in which an infected particle performs a simple symmetric random walk with jump 
rate Db, and recuperates at a given rate, while a healthy particle performs a simple 
symmetric random walk with jump rate Da- (Healthy particles are infected on contact 
with infected ones.) The ARW model corresponds to Da = 0. The generalized CP (in 
the case Da = Db) was proposed in the late 1970's by Spitzer, and later was studied 
in detail in [32l |33l |3ll |35] . The diffusive epidemic process has also been studied via 
renormalization group and numerical simulation [361 EH SSI [201 (ZD [22l [271 [HI ■ A 
general conclusion from these studies is that there are three distinct regimes of critical 
behavior, for Da < Db, Da = Db and Da > Db- It is not yet clear whether the ARW 
model falls in the first regime, or, alternatively, that Da = marks a special case. 

The ARW with symmetric hopping is closely related to the conserved lattice gas model 
(CLG) [Sni [381 [Ml [iQl [H], the principal difference being that in the CLG a site can 
be occupied by at most one particle, while active particles are those having at least 
one occupied neighbor. The CLG and the conserved stochastic sandpile share the same 
essential features, namely, a continuous phase transition between an active and an 
absorbing state, conservation of particles, and coupling between the order parameter 
(activity) and particle density, with the particle configuration frozen in regions devoid of 
activity. There is evidence [48l [TBI HOj that the CLG and conserved stochastic sandpile 
exhibit the same critical behavior, and we should expect the same to apply to the 
symmetric ARW model. 

Numerical analysis and some general theoretical arguments suggest that the ARW model 
exhibits a phase transition in the parameters A and and that there should be two 
distinct regimes: 

i) Low particle density. There is a phase transition in A in this case, namely if A is 



3 



large enough, then system locally fixates, i.e. for any finite volume A there is almost 
surely a finite time t\ such that after this time there are no A particles within A. If A 
is small enough there is no fixation, and we expect that there is a limiting density of 
active particles in the long-time limit. 

ii) High particle density. In this case there is no phase transition. For any A > 0, the 
system does not fixate. 

The balance of this paper is organized as follows. In Section [2] we study ARW2 with 
totally asymmetric walks via mean-field theory. Section [3] contains a detailed study 
of the critical behavior for totally asymmetric walks, including power laws and scaling 
relations, via Monte Carlo simulations, and ends mentioning the symmetric case, which 
belongs to another universality class. In Section[4]we quote the few known mathematical 
results and discuss a collection of open problems. 

2 Mean-field theory 

In this section we develop a mean-field theory (MFT) for the ARW2 model defined in 
Section [TJ As is usual in this type of approach, we treat the state of each site as statis- 
tically independent. Although the discussion is formulated for the totally asymmetric 
case, this 'one-site' approximation in fact yields the same predictions for the symmetric 
version. 

For n > 1, let Pn{t) be the fraction of sites having exactly n 74-particles. We denote the 
fraction of sites occupied by an S-particle by Pi{t), and the fraction of vacant sites by 
Poit)- Normalization implies, 

oo 

P'l + = 1 

n=0 

while the particle density is 

oo 
n=0 

Let 

oo 

P* = '^npn, 

n=2 

so that the density of A-particles is pa = p* + pi = C ~ p'l ■ 

We now obtain the equations of motion for the p„, starting with n = 0. The rate of 
transitions into n = is ppi, that is, to enter the state a site must have a single 
^-particle, which leaves at rate p. The rate of transitions out of state is: 

oo 

^ npn,o + PPl,0 
n=2 

where Pn,m is the joint probability for a pair of nearest-neighbor sites j and j + 1 to 
harbor n and m A-particles, respectively. (pi,o is the joint probability for a site to be 
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empty, and its neighbor on the left occupied by a single A-particle.) The reason is that 
to exit state zero, a site must be in that state and have a nearest neighbor on the left 
with one or more particles capable of jumping onto it. 

The mean-field approximation consists in factoring all joint probabilities: Pn,m — ^ PnPm- 
Combining the rates for transitions into and out of state zero, and applying the mean- 
field factorization, we obtain: 

dpo 

— =pp^- ppQ^ 

where 

p = p* +ppi. 

Proceeding in the same manner we obtain equations of motion for the other one-site 
probabilities: 

= p{Po - Pi) + 2p2 - Pi 

at 

dp'i ~ , 

^ = qpi-pp, 

dp2 I 

= P{Pl +Pl- P2) - 2p2 + 3p3 

and for n > 3, 

^ = p{pn-l - Pn) - npn + (u + l)p„+l. (l) 

These equations conserve normalization and the particle density C- For C < 1 there is 
an inactive solution, p[ = C, po = 1 — C- 

We seek an active stationary solution by introducing a 'quasi-Poisson' ansatz, 

A" 

Pn = A. — - for n > 3. 

ni 

Substituting this hypothesis in Eq. ([1]) (with n > 4) one finds X = p. Using the equations 
for n = 3, 2, 1, and 0, we obtain p2 = Ap'^/2, pi + p'l = Ap, pi = Af? /{q + p), and 
Pq =ppi/p- Normalization then implies 



A 



q + p 

With the stationary distribution in hand, we may write 



peP 



C = p*+p,+p[= A{p)peP = ' (2) 

eP ppI {q^ p)-\ 

A plot of the active particle density Pa = C ~ p'l versus is shown in Figure [H for 
p=l/2. 
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Figure 1: Comparison of simulation (squares) and mean- field theory (solid line) for the 
stationary activity density Ca at p = 0.5. 

To locate the critical point Co we evaluate the limit of the r.h.s. of Eq. ([2]) as p — > 0, 
yielding 

Cc = l-P. 

For densities smaller than (c, P = and C = Pv i-^-' particles eventually go to sleep. 
The stationary activity density grows oc C ~ Cc for C > Cc- Numerical integration of 
the mean-field equations shows that the solution indeed converges to the stationary one 
found above. For an initial Poisson distribution, or one in which a fraction C < 1 of 
sites are singly occupied and the rest vacant, the approach to the stationary state is 
monotonic, and exponentially rapid away from the critical point. At the critical point 
the activity density decays algebraically, pa ~ t~^. These are the usual characteristics 
of a mean- field theory for an absorbing-state phase transition [44] . 

The mean-field analysis is readily extended to the two-site approximation, in which 
the dynamical variables are the joint probabilities pij for a pair of neighboring sites to 
have occupations i and j [44] . (In this case three-site probabilities are approximated 
so: Pi,j,k — Pi,jPj,k/Pj-) The pair approximation again yields Cc = 1 — P- The stationary 
density of active sites is slightly less than in the simple MF approximation (a reduction 
of about 6%, near the critical point), and the relaxation at the critical point again 
follows Pa ~ t~^. Finally, we note that for model ARWl (jump rate unity for all A 
particles, sleep rate A for isolated A particles), the mean-field analysis yields the critical 
density 

Cc = A/(l + A), 

in agreement with simulations and rigorous results (Theorem [6]) . 




6 



3 Simulation results 



We performed Monte Carlo simulations of the ARW2 model on rings of L = 100 to 8000 
sites. Stationary and time-dependent properties were determined from averages over 10^ 
- 10^ independent realizations of the process, starting from an initial configuration in 
which all particles are in state A. In simulations, we select an A-particle (a list of such 
particles is maintained), and if it is not isolated, it jumps to the right. If the selected 
particle is isolated then it goes to sleep with probability q, and jumps forward with 
probability p = 1 — q. The time increment associated with each event is At = 1/n, 
with n the number of A particles just before the event. (For C < 1, a finite system 
must eventually become trapped in an absorbing configuration. In practice, however, 
the lifetime of the quasi-stationary metastable state observed in simulations is very 
long for p > pc- It appears, moreover, that the quasi-stationary properties observed in 
simulations converge to a well-defined limit as the system size L — > oo, that is, to the 
true stationary properties of the infinite system.) 

3.1 Phase diagram 

We studied the stationary density of ^-particles pa, the moment ratio m = {Pa)/Pai 
the survival probability Ps{t), that is, the probability that not all walkers are asleep. 
The variation of the stationary activity density with walker density ( (for fixed sleeping 
probability q = 1/2) is shown in Figured) Note that the data represent extrapolations 
to the infinite-size limit based on results for systems of size 100, 200,. ..,3200. The 
simulation result is very close to, although systematically smaller than, the mean-field 
prediction. (The small but nonzero difference between simulation and MFT cannot 
be attributed to a finite size effect.) The pair approximation is in somewhat better 
agreement with simulation. For example, the stationary activity density near the critical 
follows 

Pa = B{C - Cc). 

For p = 1/2, for example, simple MFT yields an amplitude of -B = 1.70; the amplitude 
in the pair approximation is 1.60, while simulation yields B = 1.46. 

Of particular interest is the location of the phase boundary. For C ^ 1 (with p fixed 
at 1/2), the data for pa fall very nearly on a straight line that intercepts pa = at 
C = 0.5, as predicted by MFT. For any finite system size the quasi-stationary activity 
density at (" = 1/2 is nonzero, but pa approaches zero with increasing system size. 
Finite-size scaling theory [471 [T8| [19] predicts a power-law dependence of the stationary 
order parameter pa on system size along the critical line: 

Pa{pc,C,L)^L-f'/''^ (3) 

where (3 and are the critical exponents associated, respectively, with the order pa- 
rameter and the correlation length [44]. (Away from the critical line pa converges 
exponentially to its stationary value as L — > oo.) A similar picture holds for the lifetime 
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t{p, L) defined in terms of tlie survival probability Ps{t). In a finite system, for C < 1) 
we expect Pg ~ exp(— t/r). Along the critical line, finite-size scaling theory predicts 

T{p,,C,L)r.L''\\/'^^ (4) 

with the critical exponent associated with the correlation time. Rather than attempt 
a systematic justification of these scaling ideas here, we simply note that the behaviors 
implied by Eqs. ([3]) and ^ have been amply confirmed in studies of many absorbing- 
state transitions (as well as in equilibrium critical phenomena), including conserved 
stochastic sandpiles. 

Power-law scaling of pa and r with L provides an effective criterion for locating the 
critical point in simulations; we use it to determine pc for = 0.25, 0.5 and 0.75. Our 
results agree, to within statistical uncertainty, with the mean-field prediction pc = 1 — C- 
In all three cases we find 

/3/i/_L = 0.5 

to within a statistical uncertainty of less than 0.5%, strongly suggesting the value 1/2 
for this exponent ratio (see Figure [2|). The lifetime r (Figure [3]) can be fit to high 
precision with the expression r = C + c'L, where C and c' are constants, so that the 
ratio 

Finally, the fact that the order parameter pa is proportional to C~Cc near the transition 
implies the exponent value 

P = l. 

Using the finite-size scaling relations we then have 

= z/|| = 2. 

The moment ratio m has been found to take a well-defined value at the critical point 
of an absorbing-state phase transition [121 [28] . Consistent with this result we find m — > 
nic = 1.298(4) at all three values studied. We note that while /? = 1 is characteristic 
of mean-field-like transitions to an absorbing state, the values of z/_|_, and nic, are not 
typical of other known universality classes for absorbing-state phase transitions [451 [Ml 
144] . We suspect that the anisotropic dynamics underlies this difference. 

The model exhibits a somewhat different behavior at the end of the critical line, C = 1) 
p = 0. In this case an isolated particle goes to sleep at rate 1, i.e., it can never jump 
forward. Simulations reveal no quasi-stationary state at this point: the activity density 
decays to zero monotonically, for all system sizes {L = 100,. ..,800) investigated. The 
activity density again grows (x C — 1 near the transition, so that (3 retains its value of 
unity. 
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Figure 2: Stationary order parameter versus system size for C = 0.25 and p = 0.75. 
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Figure 3: Lifetime r versus system size as in Figure [21 
3.2 Approach to the steady state 

We studied two kinds of initial condition. In one, N = QL active walkers are inserted 
randomly and independently into the system; in the other (for C, = 1/2), only even- 



9 



numbered sites are initially occupied by an active walker. (We call these random and 
alternating initial configurations, respectively. In studies with random initial conditions 
each realization is performed using a different, independent initial configuration.) The 
two initial states lead to the same quasi-stationary properties, but the approach to the 
latter is different in the two cases. 

Consider first the evolution of pa (averaged over 10^ independent realizations) at the 
critical point C = P = 1/2, using the alternating initial configuration. Figure [3] shows 
that the evolution is nonmonotonic, as has been found for the stochastic sandpile (with 
symmetric dynamics) at its critical point [9]. The main graph of Figure S] shows p* = 
L^/'^ Pa{t) as a function of t* = t/L. (The definitions of the scaling variables p* and t* 
are motivated by the finite-size scaling results discussed above.) Under this rescaling, 
data for L = 200, 400, 800, 4000 and 8000 collapse onto a master curve. (The collapse 
is not perfect; the secondary maximum near t* = 0.8 becomes sharper as the system 
size is increased.) 
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FIG. 4 

Figure 4: Main graph: scaled activity density p* = L^/'^pa{t) versus scaled time t* = 
t/L, for C = p = 1/2 with an alternating initial configuration. Data for system sizes 
L = 200, 400, 800, 4000 and 8000 are superposed (sharper maximum corresponds to 
larger size). Inset: activity density versus time on log scales, L=4000. 

The inset of Figure H] shows the overall relaxation to the quasi-stationary state. The 
initial decay appears to follow a power law pa '-^ t'^ with 6 = 0.50. At absorbing state 
phase transitions one expects the scaling relation 5 = which is indeed verified if 

we insert the values /? = 1 and = 2 found above. The initial growth of the moment 
ratio is expected to follow m — 1 ~ t^/^, with the dynamic exponent z equal to the ratio 
[54]. In fact we find m — 1 oc t, consistent with the exponent ratio found above. 
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For random initial conditions the general picture is similar, although there are some 
differences in detail. The relaxation is again nonmonotonic, with a collapse of data for 
various system sizes using the scaling variables p* and t* , but the secondary maximum 
(which falls near t* = 0.7), is smooth, rather than cusp-like as for the alternating initial 
condition. The initial decay again appears to follow a power law, with 5 ~ 0.51. The 
initial growth in the moment ratio follows m — 1 ~ t^/^, but with 1/z = 1.025 rather 
than the expected value of unity. The slightly larger apparent exponents observed with 
random initial conditions may reflect corrections to scaling due to relaxation of long- 
wavelength modes present in the initial distribution (and which are strictly excluded in 
the alternating case). 

Some understanding of the relaxation may be gleaned from the spatial distribution 
of the particles. Figure [5] shows the spatio-temporal evolution of a typical realization 
at the critical point, C = p = 1/2, for random initial conditions. We see that after 
an initial transient, all of the active particles are confined to a relatively narrow band. 
During coalescence into a single band, the activity density decreases rapidly; coalescence 
appears to be irreversible. Studies of larger rings confirm these observations. The time 
for the activity to become confined to a single band grows with system size, but is 
typically smaller than L/2. (We have not determined if the time grows linearly with 
L or more slowly.) The active region propagates through the system at a steady rate, 
expanding or contracting due to intrinsic fluctuations and to the varying density of 
sleeping particles it encounters as it moves. The boundaries of the active region move 
at a speed somewhat greater than unity: in large systems the speed is 1.18 - 1.20 
(sites per unit time). Repeated encounters of the active band with regions rich in 
sleeping particles may be connected with the revivals observed in the activity density 
(Figure H]). We observe the coalescence into a single band of activity, and the same 
speed of propagation, in studies with the alternating initial configuration. (Naturally 
the initial transient is different in the two cases.) In the infinite system, activity cannot 
be confined to a narrow band, but we should expect, on the basis of the foregoing 
observations, a steady coarsening of the activity pattern. 

We have also studied the distribution of first passage times tq to the origin of the ARW2 
model on the line (i.e., in a system without periodic boundaries); tq is defined as the time 
at which a particle first jumps from site -1 to the origin. To study its distribution we 
simulate the system on the lattice extending from x = —L to x = 0. We determine the 
probability density p{tq) up to a certain maximum time, by studying a series of lattice 
sizes L, until p{tq) stabilizes. Figure [6] shows the density obtained for a lattice size of 
10^ sites, for a system at the critical point, C, = p = 1/2. The same result is obtained 
for 10^ sites, within the uncertainty. The data can be fit with a power-law, p{tq) ~ Tq"") 
with a = 1.50(1). Thus the mean first passage time to the origin diverges at the critical 
point. It is curious that an interacting particle system with totally asymmetric jumps 
shares the same exponent as that of a symmetric random walk. For an unbiased random 
walk the first return is the smallest time when it is positive, whereas tq is determined 
by the smallest lattice interval [— x, 0] where the presence of initial particles exceeds the 
gaps of activity. This subtle heuristics explains why the tail of tq decays with exponent 
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Figure 5: Typical evolution of a system at the critical point C = P = 1/2, with a random 
initial configuration, L = 160. In each horizontal sweep, the height of the line represents 
the number of particles at site x, with a sleeping particle corresponding to height -1. 
The graph at the left shows the number n of active particles versus time. 

1/2, though it is not so clear why its density is so smooth. 

3.3 The symmetric case 

We performed a series of studies of the symmetric ARW2 model at density C = 1/2. 
In contrast to the asymmetric case, the critical value of the hopping probability, pc, is 
considerably higher than its mean-field value, Pc = min{^, 1}. In terms of the ARWl, 
it means that 

A 

^'^^TTa' 
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Figure 6: Distribution density of the time of the first jump to the origin for a system 
at the critical point C = P = 1/2, with an alternating initial configuration, and L = 10^ 
sites. The blue line has slope —1.50. 

and there are partial mathematical results in this direction (Theorem [5]) . Studies using 
rings of up to 12800 sites yield pc = 0.87835(2). (Our criterion for criticality is that 
the stationary activity density follow a power law, pa ~ L~f^/'^^.) These studies yield 
the exponent ratio (3/v±_ = 0.23(1), very different from the result of 0.5 found in the 
asymmetric case. A study of the growth of m — 1 at the critical point, in a system 
of 12800 sites, yields the dynamic exponent z = 1.51(1). The stationary value of the 
moment ratio m at the critical point is rric = 1.15(1). The results for /9/i^_l, z and 
rric are all quite far from the corresponding values in the asymmetric case. They are, 
on the other hand, rather close to those found for a conserved stochastic sandpile |10j : 
/3/i/_L = 0.217(6), z = 1.50(4) and = 1.14(1). 

These results support the assertion (based on considerations of symmetry) that the 
symmetric ARW model falls in the conserved stochastic sandpile universality class. Since 
scaling properties of sandpile models are rather subtle, we defer a full characterization 
of the symmetric model to future work. It is nevertheless clear that the symmetric 
and asymmetric ARW exhibit very different critical behavior. On a qualitative level 
the difference is quite dramatic if we compare the evolution of the asymmetric model 
(Figure [5|) with that of the symmetric model at its critical point (Figure [7|). In the latter 
case there is no tendency for the activity to become confined irreversibly to a narrow 
band; active regions are seen to branch as well as coalesce. 

As noted above, the model studied by Jain [26] corresponds to setting p = in the 
symmetric ARW. The properties demonstrated in p6] are in fact very different than 
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Figure 7: A typical realization of the symmetric ARW model near the critical point. 
Time increases downward, with white, light blue, dark blue and black points represent- 
ing empty sites, sites with an inactive particle, sites with an active particle, and sites 
with > 2 particles, respectively. 

those obtained here, for p > 0. In the former case the active phase, which exists for 
C > Cc = 1) has Pa oc C~Cc, so that the critical exponent /3 = 1, and the stationary prob- 
ability distribution in the active phase is uniform on the set of allowed configurations 
(i.e., those in which no site is empty). Since the active phase has a product measure, 
correlation functions are identically zero and the critical exponent i'± is undefined. The 
difference between the scaling properties observed for p > 0, and those found in [26] 
may be understood by noting that in our case an isolated particle, while active, may 
evade becoming immobile by jumping to an occupied site. In this process it may re- 
activate a sleeping particle. Thus the number of active particles fluctuates while total 
particle number is conserved, a hallmark of models in the conserved stochastic sandpile 
universality class. When p = 0, by contrast, the number of mobile particles is fixed at 
— L in the stationary state. 

It is also worth noting that although the asymmetric ARW (with p > 0) shares the 
critical exponent value (3 = 1 with the model studied in [26], it is different since the 
exponent is well defined, and the model satisfies finite-size scaling. Again, the 
asymmetric ARW features a fluctuating number of active particles when p > 0. 
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4 Rigorous results, conjectures and open problems 



In this section we will summarize the few existing rigorous results concerning the ARW 
model, discuss some conjectures and present several open problems, whose understand- 
ing may shed some light on the long-time behavior of the system. 

A rigorous understanding of this model is still in its embryonic stage, and some of the 
open questions appear to be quite difficult and mathematically challenging. 

4.1 General case 

We start with the first basic fact, proved in |49j using the Diaconis- Fulton [TJ [T7] rep- 
resentation of the model. The representation provides an Abelian property for the dy- 
namics of the system with finitely many particles, and - what is particularly important 
- provides monotonicity for the occupation times in as well as in A. 

Theorem 1 ([IH]). For d ^ 1 and any translation- invariant random walk and A > 0, 
there exists = Cc(A) € [0, oo], such that if the initial distribution is i.i.d Poisson with 
density then 



Moreover, (^c is non-decreasing in A. 

For fixed A the value of Cc(A) is not known, however some theoretical arguments suggest, 
and numerical simulations support, that the following holds: 

Conjecture 1. For any dimension, any random walk, and any A > 0, 



4.1.1 Supercritical regime 

Using Peierls type argument one can show that Cc(A) < +oo: 

Theorem 2 ([M]). Consider simple symmetric random walks on ^ 1. There 

exists Co < oo such that Cc(A) < Co for all A. 

Recently E. Shellef improved this estimate: 
Theorem 3 ([53]). Under the same hypotheses. 





C<Cc 
OCc- 



< Cc(A) < 1. 




Cc(A) ^ 1. 



(5) 



Another approach to prove ^ is to show mass conservation for this model: 
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Theorem 4 ([T]). For i.i.d. initial conditions, simple symmetric random walks, if there 
is local fixation, then each particle jumps finitely many times. By the mass transport 
principle, this implies that whenever there is fixation the density of the limiting state 
in the same as the density of the initial state, in particular it implies ([5]). 

However, the following problem remains open: 
Problem 2. Show strict inequality in 

Another interesting question about the supercritical regime is the following: 

Problem 3. For each ^ > i^c, show that there is a unique non-trivial invariant distri- 
bution, ergodic with respect to spatial translation, whose particle density is 

4.1.2 Subcritical behavior 

On the other hand, it is rather easy to get convinced that Q ought to be strictly positive. 
Despite the recent progress in the one-dimensional case (see [49]), there are no results 
in higher dimensions: 

Theorem 5 ([M]). For d = 1, bounded range random walks and any A > we have 
that < Cc ^ 1- For nearest-neighbor walks we have: 

^ ^ Cc ^ 1. 



1 + A 



The proof of the above theorem again relies on the Diaconis-Fulton representation of 
the dynamics, in particular it uses the Abelian property and certain monotonicity for 
finite particle systems. Conceptually all the ingredients of the proof can be used in 
higher dimensions, with the exception of the last estimate, which in dimensions ^ 2, 
if repeated straight forwardly, boil down to the need for refined bounds on growth 
intensity of a Diffusion Limited Aggregation type growth model, and the analog of the 
one-dimensional argument produces unsatisfactory estimates. Thus, we have 

Problem 4. Show that Cc > for d ^ 2. 

Though Theorem [5] establishes the fact that (^c > in one-dimension for a rather broad 
class of walks, it does not give a satisfactory description of the final absorbing state, 
nor correct estimates for the fixation time. We may therefore state several important 
questions. 

Problem 5. Describe the distribution of the final configuration after fixation (in any 
dimension, including d = 1). Numerical analysis suggests that it depends on the initial 
distribution. It seems however to be less sensitive when ^ ~ (^c, when the final state is 
prominently diff'erent from a Bernoulli. 
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Problem 6. Establish the rate of fixation, i.e., the asymptotic behavior of the proba- 
bility that there is an active particle at the origin after time t, for t large enough. We 
believe that if C is small enough, the decay should be exponential. However when ap- 
proaches (c we do not exclude the possibility that decay may become slower (stretched 
exponential, or even algebraic). 

The following problems constitute possible intermediate steps to understand the behav- 
ior of the system in the subcritical regime. 

Problem 7. Consider an infinite volume system with only k < +00 particles, which 
are all initially located at the origin and in state A. Let Tfc denote the (a.s. finite) time 
when the system fixates, i.e., when all particles become inactive. For any fc > find an 
asymptotic bound for P{Tk > n) as n ^ +00. Prove the Large Deviation Principle for 

n- 

Problem 8. As before, consider a finite system of k particles and denote by Lk{s) and 
Rk{s) positions of the leftmost and rightmost particles in the system, and by L^{s) 
and R^{s) position of the leftmost and rightmost ^-particles in the system - which 
do not necessarily coincide with Lk{s) and Rk{s), but for all times satisfy inequalities 
L^is) > Lfc(s) and R^{s) < Rk{s). 

How does R^{s) — L^{s) behave during the time interval [0,Tfc]? 

What can we say about the distribution of Rk{Tk) — Lk{Tk), the diameter of the con- 
figuration in the final state, and about its displacement with respect to the origin 
\0-Lk{Tk)\? 

4.1.3 At criticality 

We have good reason to believe that at the critical density the system does not fixate in 
any dimension (see comments in the next subsection). However time intervals between 
successive visits of active particles to the origin will diverge to infinity. It is not clear 
that here we are observing an aging phenomenon in one of accepted senses. 

4.2 Totally asymmetric dynamics in one dimension 

The description of the totally asymmetric walk in one dimension is relatively well un- 
derstood. Let xo{t) denote the position of a tagged particle starting at the origin. 

Theorem 6 ([25]). For d = 1, and the totally asymmetric walk. 
If C < Cc, then 

P[r/f (0) > for some s > t] ^ cie^^'^i^ 
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If C > Cc, then 

lim ^ = v,>0. (6) 

t — ^oo t 

If = ^c, and initially all particles in the system are A-particles, then the system does 
not fixate (!). However, 

hm ^ = 0. (7) 

t^oo t 

The limit ([6|) implies that there is a current and that the tagged particle has an asymp- 
totic velocity. This is related to Problem [3l as the limiting speed should be the density 
of active particles in the (unique) non-trivial ergodic invariant distribution with total 
density of particles (. 

The limit ([7]) tells us that there is no current in the system. This should be related to 
the absence of a non-trivial ergodic invariant distribution with total density Co i-e., the 
phase transition is not of first order. 

Moreover, it motivates the following problem. Consider a system with critical density, 
and denote by Tj(0) the time spent between arrivals of the (i — l)-th and i-th ^d-particle 
at the origin. 

Problem 9. Characterize the behavior of rj(0) when i — > +oo. 

We believe that the system exhibits behavior reminiscent of so-called aging for disor- 
dered systems. 

Problem 10. Explain the origin of the 'spikes' in the time-dependent density of active 
particles in Figure H]) . 
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